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We present measurements of the linear polarization amplitudes and the strong relative phases 
that describe the flavor- untagged decays B^ — > J/ipK*° and B 9 — > J/ip</> in the transversity basis. 
We also measure the mean lifetime r s of the B 9 mass eigenstates and the lifetime ratio r s /r d . The 
analyses are based on approximately 2.8 fb _1 of data recorded with the DO detector. From our 
measurements of the angular parameters we conclude that there is no evidence for a deviation from 
flavor SU(3) symmetry for these decays and that the factorization assumption is not valid for the 
B° d -> J/ipK* 9 decay. 
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B mesons are fertile ground to study CP violation and 
search for evidence of new physics. There are elements, in 
addition to CP violation, involved in the theoretical de- 
scription of B meson decays, such as flavor SU(3) symme- 
try, factorization and final-state strong interactions. To 
understand the role CP violation plays in these decays, 
it is essential to understand and isolate the effect of each 
of these elements in the B meson decays. 

Factorization states that the decay amplitude of B 
meson decays can be expressed as the product of two 
single current matrix elements [l[ and this implies that 
the relative strong phases are (mod it) [2|. A different 
measured value for the strong phases would indicate the 
presence of final-state strong interactions. The B® meson 
can be formed by replacing the s quark with the d quark 
in the B® meson. From flavor SU(3) symmetry applied 
to the B®-B® system one expects that the theoretical 
description is similar; in particular the B® —> J/ipK*° 
and Bg — > J/ipfj) J31 decays, can be described in the 
transversity basis [2| by the linear polarization ampli- 
tudes, Aq, An, and A±, and the relative strong phases Si 
and 62- Flavor SU(3) symmetry requires that the ampli- 
tudes and phases characterizing these decays should have 
the same values. 

Other observables of these decays are the lifetimes of 
both mesons, which allow us to compare with theoretical 
predictions of the lifetime ratio. Phenomenological mod- 
els predict differences of about 1% 0, [1| between the B^ 
and B® lifetimes. Previous B meson lifetime measure- 
ments [6j are consistent with these predictions. 

In this Letter we report the measurements of the pa- 
rameters that describe the time-dependent angular dis- 
tributions of the decays B Q d -> J/tpK*° and B° -> J/<ip<f> 
in the transversity basis, where the initial B meson fla- 
vor is not determined ("untagged"). We study the B® 
and B° a mesons to verify the validity of the factorization 
assumption Q and to check if flavor SU(3) symmetry Q 
holds for these decays. We also report the lifetime ra- 
tio f s /Td for these mesons and the width difference Ar s 
between the light and heavy B® mass eigenstates. The 
analyses were performed using data collected with the DO 
detector 0] in Run II of the Fermilab Tevatron Collider 
during 2003 — 2007 with an integrated luminosity of ap- 
proximately 2.8 fb _1 of pp collisions at a center-of-mass 
energy of 1.96 TeV. In contrast with the flavor-tagged 
analysis reported in Ref. Q, in this Letter we report a 
simultaneous analysis of both the B® and B® meson de- 
cays, carried out in such a way that a straightforward 
comparison between their angular and lifetime parame- 
ters can be performed. 

We use the B° s -► J/V>&J/V> -> 4> -> K+R- 

selection described in Ref. jg. The decay B° -> J/ipK*°, 
J /if) — > K*° — > K ± n T is reconstructed using 



similar selection criteria and algorithms as the B® chan- 
nel because they have the same four-track topology in 
the final state. The differences are the requirement that 
the transverse momentum of the pion be greater than 
0.7 GeV/c, the invariant mass for the ( J/ip, K*°{892)) 
pair be in the range 4.93 — 5.61 GeV/c 2 , and the se- 
lection of the if* (892) candidates by demanding the 
two-particle invariant mass between 850 MeV/c 2 and 
930 MeV/c 2 . Due to lack of charged particle identifi- 
cation, we assign the mass of the pion and kaon to the 
latter two tracks and use the combination with invariant 
mass closest to the K*° mass. 

The proper decay length (PDL), defined as in Refs. 0, 
ll| , for a given B® or B° a candidate is determined by mea- 
suring the distance traveled by each 6-hadron candidate 
in a plane transverse to the beam direction, and then 
applying a Lorentz boost correction. In the B® and B® 
final selection, we require a PDL uncertainty of less than 
60 fim. We find 334199 and 41691 candidates that pass 
the £?2 and B® selection criteria, respectively (see Fig.Q]). 

We denote the set of the angular variables defined in 
the transversity basis, where the decays B® — ► J /ipK*° 
and B® — > J/tpcj) are studied, as u> = {cp, cos 9, cos tp}. 
The description of these decays in this basis gives us ac- 
cess to the three linear polarization amplitudes at pro- 
duction time, t = 0, |A (0)|, I -Ay (0)1, and |Al(0)|, satisfy- 
ing I^oP + IAhP + IAJ 2 = 1 42]; and the CP-conserving 
strong phases Si = argL4jMx], and S 2 = arg[AgAx]- 
Since only the relative phases of the amplitudes can enter 
physics observables, we are free to fix the phase of one of 
them, and we choose to fix S = arg(Ao) = 0. 

According to the standard model, CP-violation effects 
in the B® system are very small [13(. In this analysis, 
we assume CP conservation and express the differential 
decay rate for the untagged decay B° — > J/ip4> as Q: 

d*V/(dwdt) oc e- r ^ [\A \ 2 fi(ij) +Re(A*A ll )f 5 (ij) 

+ |A,||V a (uO]+e- rH Vx| 2 /3(w)> (1) 

where ^l(H) = 1 I t l{h) is the inverse of the lifetime corre- 
sponding to the light (heavy) mass eigenstate. The mea- 
sured parameters, the width difference Ar s = Tl — r# 
and the mean lifetime f s = 1/T = 2/ (r& + Th), are 
given in terms of these inverse lifetimes. The angular 
functions /i(u>) are defined in Ref. Q. In this decay, we 
have access to the phase S\\ — arg(Aoj4ii), which is related 
to S± and S2 by 5\\ = S2 — Si. 

In the B® system, there is evidence of interference be- 
tween the P- and S- wave Kir amplitudes [14j |. which is 
taken into account in this analysis. The differential de- 



cay rate for the untagged decay B® 
by@,[l4|: 



J/ipK*° is given 



d 4 V/(du>dt) oc e - rdt {cos 2 A[|A | 2 /i(w) + |^||| 2 / 2 (u;) 
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\A ± \ 2 Mu)-CM^±)M") 
Re(4SA||)/ 5 (w) + C Im(ASA x )/6(w)] 
sin 2 A ■ /7(w) 



•sin2A [/ 8 (w)cos (S {] - S s ) \A\\\ 



+ f g (u)sm(6±-8 a )\Ai 
+ /i (w)coB<y»-|^o|]}, 



(2) 



where r,j = I/tj is the inverse of the lifetime, £ = 
+1(C = —1) for K + (K~); A, <5 S , and are defined in 

Refs. @,[I3. For the Ar d is expected to be zero [iljj . 

An unbinned likelihood fit is performed to extract all 
the and i?" parameters. For the j'th B meson can- 
didate, the inputs for the fit are the mass m,j, PDL 
PDL uncertainty er ct . , and the angular variables 
The likelihood function £ for the untagged decays 
-> J/i[)K*° and B" -> is defined by 



£ = nti + (i - /. 



(3) 



where iV is the total number of selected events and f s is 
the fraction of signal events in the sample, a free param- 
eter in the fit. 

T s is the product of the signal probability distribution 
functions (PDF) of mass, PDL, and transversity angles, 
and the angular acceptances, which are determined via 
Monte Carlo simulations. The mass and PDL signal dis- 
tributions are modeled for both decays in the same way. 
The mass distribution is modeled by a Gaussian func- 
tion with free mean and width. The PDL distribution 
is described [T3| by the convolution of an exponential, 
whose decay constant is one of the fit parameters with 
a resolution function represented by two weighted Gaus- 
sian functions centered at zero. The widths Si<x ct . of each 
Gaussian with scale factors Sj (i = 1,2) are free param- 
eters in the fit to allow for a possible misestimate of the 
PDL uncertainty. The transversity angular distributions 
are modeled by the corresponding normalized Eqs. ([JJ) 
and |(3J) . The contribution where the mass of the K and it 
are misassigned in our data is estimated by using Monte 
Carlo studies to be about 13% and is taken into account. 

!Ft is the product of the background PDF of the same 
variables and the angular acceptance as in the signal. We 
separate the background contributions into two types. 
The prompt background accounts for directly produced 
J ftp mesons combined with random tracks. Non-prompt 
background is due to J / ip mesons produced by a b hadron 
decay combined with tracks that come from either a 
multibody decay of the same b hadron or from hadroniza- 
tion. The mass distribution for the background is mod- 
eled by two independent normalized negative-slope ex- 
ponentials, one for the prompt and one for the non- 
prompt contributions. The PDL distribution for the 
prompt background is parameterized by the resolution 
function described above. The PDL distribution for the 



non-prompt background is modeled by a sum of two ex- 
ponential components for positive ct and one for negative 
ct that account for a mix of heavy flavor meson decays 
and their possible misreconstruction. The angular distri- 
butions for the background components are modeled by 
a shape similar to that of the signal, but with an inde- 
pendent set of amplitudes and phases. 

The results of our measurements are summarized in 
Table |TJ Figures Q] and [5] show the mass and the PDL 
distributions for the B^ and B® candidates, respectively, 
with the projected results of the fits. The parameters 
with the strongest correlations are the linear amplitudes 
for the Bj , and the width difference and the mean lifetime 
for the B°. 



TABLE I: Summary of measurements for the decays S° — > 
J/tpK*° and B® — > J/t/j<j>. The uncertainties are only statis- 
tical. 



Parameter 


B° d 


B° 3 


Units 


\Ao\ 2 


0.587 ±0.011 


0.555 ±0.027 




IV 


0.230 ±0.013 


0.244 ± 0.032 




Si 


-0.38 ±0.06 




rad 


S 2 


3.21 ±0.06 




rad 


h 




9 79+112 
/5 - ' ^-0.27 


rad 


T 


1.414 ±0.018 


1.487 ± 0.060 


ps 


AT S 




085+ 072 

U.U8O_ 07g 


ps" 1 


N 3i g 


11195 ±167 


1926 ± 62 





Table |TT] summarizes the systematic uncertainties in 
our measurements for B® and B® decays. To study the 
systematic uncertainty due to the model for the mass dis- 
tributions, we vary the shapes of the mass distributions 
for background by using two normalized first-order poly- 
nomials instead of the nominal two negative exponentials. 
We estimate the systematic uncertainty due to the res- 
olution on the PDL by using one Gaussian function for 
the resolution model. The fitting code is tested for the 
presence of biases by generating 1300 pseudo-experiments 
for £>J] and 1000 for B° s , each with the same statistics as 
our data samples. We generated the events following the 
PDL, mass, and transversity angular distributions de- 
scribed above. The differences between the input and 
output values are quoted as the systematic uncertainty 
due to the fitting. The systematic uncertainty for 5\\ re- 
ported for this source is due to an intrinsic ambiguity 
for this parameter in Eq. (TT]). The pseudo-experiments 
produced also cover the other solution for Su . The contri- 
bution from the detector alignment uncertainty is taken 
ll| . Other potential sources of systematic 



from Ref. 

uncertainties have been investigated and found to give 
negligible variations in the measured parameters. The 
systematic uncertainties for the ratio f s /Y<2 are obtained 
by finding the ratio of the lifetimes for each systematic 
variation on Table [TT] and taking the difference between 
this value and the nominal ratio. 
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TABLE II: Summary of systematic uncertainties in the measurement of angular and lifetime parameters. The total uncertainties 
are given combining individual uncertainties in quadrature. 





B° s 




Source 






Si (rad) 


^2 (rad) 


Td (PS) 


\A \ 2 


\ A \\\ 2 


$11 (rad) 


Ar s (ps- 1 ) 


T 3 (PS) 


T s /r d 


Mass background 




0.024 


0.09 


0.05 


0.030 


0.004 


0.002 


0.02 




0.021 


0.009 


PDL resolution 


0.013 


0.008 


0.02 


0.03 


0.013 


0.005 


0.003 






0.016 


0.012 


Fitting code 


0.001 








0.004 


0.004 


0.014 


0.26 


0.001 


0.008 


0.003 


Alignment 










0.007 










0.007 




Total 


0.013 


0.025 


0.09 


0.06 


0.034 


0.006 


0.014 


0.26 


0.001 


0.028 


0.015 



DO, 2.8 fly 1 




&> 
a 

S3000 

B 

> 

w 2000: 



1000: 



Total fit 

Prompt background 

Non prompt background 



5.1 5.2 5.3 



4) 

s 

u 

v 
a 

a 

s 

> 



800 
700 
600 
500 
400 
300 
200 
100 




5.4 5.5 5.6 

Mass (GeV/c 2 ) 

DO, 2.8 fly 1 




Total Fit 

Prompt background 

Non prompt background 



5.1 5.2 5.3 5.4 



5.5 5.6 5.7 

Mass (GeV/c 2 ) 



© 

=>io" 

M 

0) 

Sio 3 

h 

10 

i 



L(a) A 


DO, 2.8 fb 1 


j I • 


Total fit 

Signal 

t Background 

-Tv*. 


r / 




f 


fliffe^ 


1 1 





-0.1-0.05 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

ct (cm) 



DO, 2.8 fly 1 



Total fit 




Background 



-0.1-0.05 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

ct (cm) 



FIG. 1: Invariant mass distribution for selected (a) and 
(b) B® candidate events. The points with error bars represent 
the data, and the curves represent the fit projections for the 
total and the background components. 



FIG. 2: PDL distribution for selected (a) B% and (b) B° 
candidate events. The points with error bars represent the 
data, and the curves represent the fit projections for the total, 
signal, and background components. 



In conclusion, we have measured the angular and 
lifetime parameters for the time-dependent angular un- 
tagged decays B^ -> J/ipK*® and B° s -> J/ip<p, the 
lifetime ratio of both B mesons, and the width differ- 
ence Ar s for the B® meson. From the measured life- 
time parameters f s and we obtain the ratio Tg/r^ = 
1.052 ± 0.061 (stat) ± 0.015 (syst) which is consistent 
with the theoretical prediction [5j and previous mea- 
surements @. The measurement of the width difference 
Ar s = 0.085±!j;^8 (stat) ±0.006 (syst) ps" 1 is consistent 
with the theoretical prediction [5j, [13J and with the value 
reported in Refs. 0,0]. DO also has a measurement of 
Ar s in a flavor-tagged analysis of B® — > J/ip't 1 m Re£ Q. 

Our measurements for the linear polarization am- 



plitudes for the B®, taking into account the inter- 
ference between the Ktt S-wave and P-wave, are 
\A Q \ 2 = 0.587 ± 0.011 (stat) ± 0.013 (syst) and |A,|| 2 = 
0.230 ± 0.013 (stat) ± 0.025 (syst); and for B° s : \A \ 2 = 
0.555 ± 0.027 (stat) ± 0.006 (syst), and |A||| 2 = 0.244 ± 
0.032 (stat) ± 0.014 (syst) are consistent and competitive 
with those reported in the literature @, Q, EH • Our mea- 
surement of the strong phases Si and 82 indicates the 
presence of final-state interactions for the decay B® — > 
J/1PK* since 5i = -0.38±0.06 (stat) ±0.09 (syst) rad 
is 3.5cr away from zero, where a is the total uncertainty. 
From the comparison of the measured amplitudes and 
strong phases [17[ for both decays we conclude that they 
are consistent with being equal for B® and B® and hence 
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there is no evidence for a deviation from flavor SU(3) 
symmetry. In our sample we find that the Kn S-wave 
intensity, as described in Ref. [14j| . is (4.0 ± 1.0)%. 
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